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ABSTRACT 

We study the problem of scalar particle production after inflation by an in- 
flaton field which is oscillating rapidly relative to the expansion of the universe. 
We use the framework of the chaotic inflation scenario with quartic and quadratic 
inflaton potentials. Particles produced are described by a quantum scalar field Xi 
which is coupled to the inflaton via linear and quadratic couplings. The particle 
production effect is studied using the standard technique of Bogolyubov transfor- 
mations. Particular attention is paid to parametric resonance phenomena which 
take place in the presence of the quickly oscillating inflaton field. We have found 
that in the region of applicability of perturbation theory the effects of parametric 
resonance are crucial, and estimates based on first order Born approximation often 
underestimate the particle production. In the case of the quartic inflaton poten- 
tial V((p) = Ay? 4 , the particle production process is very efficient for either type 
of coupling between the inflaton field and the scalar field x even for small values 
of coupling constants. The reheating temperature of the universe in this case is 
[A log (1/A)] times larger than the corresponding estimates based on first order 
Born approximation. In the case of the quadratic inflaton potential the reheating 
process depends crucially on the type of coupling between the inflaton and the 
scalar field x an d on the magnitudes of the coupling constants. If the inflaton 
coupling to fermions and its linear (in inflaton field) coupling to scalar fields are 
suppressed, then, as previously discussed by Kofman, Linde and Starobinsky (see 
e.g. Ref. 13), the inflaton field will eventually decouple from the rest of the mat- 
ter, and the residual inflaton oscillations may provide the (cold) dark matter of 
the universe. In the case of the quadratic inflaton potential we obtain the lowest 
and the highest possible bounds on the effective energy density of the inflaton field 
when it freezes out. 



1. Introduction 

According to the simplest version of the inflationary scenario (see [1]) the 
universe in the past expands almost exponentially with time (such an expansion is 
called "inflation") while its energy density is dominated by the effective potential 
energy density of a special scalar field tp, called the inflaton. Sooner or later 
inflation terminates and the inflaton field starts quasiperiodic motion with slowly 
decreasing amplitude. Right after inflation the universe is empty of particles, i.e. 
"cold". Quasiperiodic evolution of the inflaton field leads to creation of particles 
of various kinds, after thermalization of which due to collisions and decays the 
universe becomes "hot" . The role of the inflaton field in the inflationary scenario is 
not necessarily played by a fundamental scalar field. It can be the expectation value 
of a scalar operator constructed of fields of other type. It can also be the effective 
scalar introduced in the inflationary models based on a higher derivative theory 

[21 

of gravity. Almost all such scenarios have an important feature: the inflationary 
stage is followed by quasiperiodic evolution of the effective inflaton field that leads 
to particle creation. 

In this work we are going to study the transition of the universe from an 
inflationary to a hot stage sketched above. The problem is not new, and there 
are many papers devoted to its study using different methods in the context of 
various inflationary scenarios (see also Refs. 1 and 9). There is nevertheless 

a drawback common to most of them - the calculations were based on ordinary 
perturbation theory. The rates of particle production by the oscillating inflaton 
field were calculated in first order Born approximation. However, as it has been 
shown in Ref. 10 this approach disregards possible parametric resonance effects 
which can enhance the rate of boson particle production. To detect such effects 
one has to go beyond standard perturbation theory. Such resonance amplification 
effects have been studied in [10] in the context of the new inflationary scenario. 
It was shown that typically the oscillating inflaton field produces many particles 
due to parametric resonance effect even in the case of extremely small couplings, 



and that estimates based on ordinary perturbation theory miss the largest effect. 
Parametric resonance effects as relevant to reheating of the universe were also 
studied in Ref. 11 where it was also demonstrated that there is no resonance 
enhancement of fermionic particle production. 

In the present paper we are going to study the problem in a systematic way. We 
will consider a model with inflation based on a power law potential for the inflaton 
field. The inflaton field (p throughout this paper will be regarded as classical. In 
Section 2 we will describe the evolution of the inflaton field after inflation without 
taking into account its couplings to other fields. Couplings of the inflaton field (p to 
fields describing fermionic spin-1/2 and bosonic scalar particles will be considered 
in the following sections. In particular, the process of production of these particles 
by the oscillating inflaton field will be studied. In all the cases we will assume the 
couplings to be sufficiently small so that perturbation theory is valid. In Section 
3 we will derive the expressions for particle production rates in first order Born 
approximation of ordinary perturbation theory. The material of Sections 2 and 3 
will closely follow a previous work [12]. In Section 4 we study general aspects of the 
parametric resonance effect, and in Sections 5 and 6 we apply the results obtained 
to the concrete cases of quartic and quadratic inflaton potentials. In Section 7 we 
present our general conclusions. 

It will be shown that in many important cases it is not legitimate to neglect 
parametric resonance effects. These effects lead to a great enhancement of the 
particle production as compared to the estimates based on ordinary perturbation 
theory. The number of particles N^ in a mode with wavenumber k that has ex- 
perienced the resonance turns out to be much larger than unity. The universe 
reheating temperature therefore can turn out to be sufficiently large even for small 
values of the couplings between inflaton and other fields. 

The opinion of the authors of Ref. 11 is that parametric resonance effects of 
scalar particle production are likely to be suppressed by the processes of scattering 
and decay of the particles produced. The reason is that particle scattering and 



decay will tend to reduce the mean occupation numbers N^ in each particular mode 
so that these numbers will grow not so fast as to be able to become large during the 
resonance period. This question is very important and deserves thorough study. 
We do not aim to touch on it in this work, but hope to consider it in future 
publications. 

In parallel to our own work, Kofman, Linde and Starobinsky have also been 
working on these issues for several years. Their results were recently summarized 
in Ref. 13. In Section 7 we will compare our results with those of Ref. 13. 

2. Inflaton field dynamics 

In this section we will derive approximate analytic expressions for the evolution 
of the inflaton field as it undergoes quasiperiodic motion in the expanding universe, 
and compute the resulting effective equation of state of the inflaton stress-energy 

Consider the chaotic inflation scenario based on the scalar (inflaton) field dy- 
namics. The Lagrangian of the model is 

1 M 2 

L = -(V V ? -V(<p)-^R, (1) 

with the inflaton field potential 

V(<p) = \v*- q \<p\ q , (2) 

where q is an arbitrary positive power. We are working in the standard system of 
units in which h = c = 1, and Mp denotes the Planck mass. 

The equations of motion for a homogeneous isotropic universe and for a homo- 
geneous scalar field look like 



(p + 3#y? + V'(tp) = , (4) 

where H = a /a is the Hubble parameter, k — 0, ±1 corresponds to different signs 
of spatial curvature, dots denote derivatives with respect to the cosmological time 
t, and prime denotes derivative with respect to tp. 

From the equation (4) it follows that 

= a ~ 3 (const - / V'(tp)a 3 dt J . (5) 

This equation describes both the rapid evolution of tp, for which (p oc a~ 3 , and the 
regime of slow rolling-down during inflation, when tp ps — (1/3H) V'(tp). 

According to the inflation scenario, during inflation the scalar field tp is rolling 
down the slope of its potential from its relatively large value towards its minimum. 
The conditions for slow rolling and for inflation are, respectively, 

\tp\^3H\tp\. (6) 

and 

\H\<^H 2 , (7) 

During inflation the term k/cl 2 in the left-hand side of Eq.(3) soon becomes 
insignificant, and the conditions (6) and (7) read, respectively, 



(aA 7 m) // «^5aA 7 m 



(V'( V )) 2 « 2 ^V\ V ). (9) 

For the scalar field potential (2) these two conditions become respectively 

.p*»^M$, (10) 



V, 2 »^M2. (11) 

Note that conditions (10) and (11) are essentially the same. 

As the scalar field evolves towards its smaller values, the conditions (10) and 
(11) cease to be valid, and a new regime for the scalar field begins, namely, the 
regime of quasiperiodic evolution with decaying amplitude. To describe this new 
cosmological period, it is convenient to rewrite the system of equations (3), (4) 
(with k = 0) as follows 

p = -3Hip 2 , (13) 

where 

P = \<P 2 + V(cp) (14) 

is the scalar field energy density. Introducing the positive value tpo(t) by the 
relation 

V(<po(t)) = Pit), (15) 

we can present the evolution of the scalar field in the form 

ip(t) = <p Q (t) cos / W(t)dt , (16) 



where W(t) is some unknown function of time. Note that the representation (16) 
for the scalar field evolution with (p$ defined by (15) is always possible. Using the 
equations (2) and (12)-(15), one can then derive the following exact expression for 
the function W(t) 




The signs "+" and "— " correspond respectively to the cases (p > and (p < 0. The 



second term in the brackets in (17) is much less than unity if 



,;-2 



*«£{*¥)' *■ (18) 

This condition is just the opposite of (10), (11). When (18) is valid, the expression 
(17) for W(t) acquires the simple approximate form 



W « ,pf^>, (19) 

and, in addition, the following condition becomes valid 

\<Po/<Po\ < W, (20) 

which allows us to regard the evolution of tp(t) as quasiperiodic with slowly decaying 
amplitude <po(t). Henceforth we assume that the condition (18) (hence also (19) 
and (20)) is satisfied. 

Our aim now will be to derive approximate evolution equations for the values 
p(t) and tpo(t) under the condition (18). First we average the equation of motion 
(13) over the time period T of quasiperiodic motion of tp. Rewriting (13) in accord 
with (14) as 

P = -6H(p-V(<p)), (21) 

and taking the average of both sides we get 

T 
^P = ~JH(p-V&))dt, (22) 

o 

where Ap(T) is the change of p over the time period T. Now take into account 
that according to (20) and to the equations (15) and (12) the values tpo(t), p(t) 
and H(t) change only insignificantly during one period of oscillation of the ^-field. 



This enables us to replace all variables except (p under the integral in (22) by their 
averaged values. We will have then, to a good approximation 

T 



Ap(T) _ 6 



^Hj(p-V(<p))dt, (23) 



o 

where H and p in (23) and below denote the corresponding averaged values. The 
integral in (23) can be evaluated as 



T 



s yft=vm> \ s ^m 



1 / , - . . -<fio —(pa 



, (p-VM)* « =^— * ^F— » -±f P . (24) 

± J r dip r dcp H ' * 

-i v 7 ^^) -i V^nd 

In (24) we took into account the relation (14), integrated in the numerator by 
parts and made use of the equation (2). In the approximation considered we 
can replace the finite difference expression on the left-hand-side of (23) by the 
derivative. Equation (23) then becomes 

p « ~^-ffp. (25) 

All the values in (25) are now to be regarded as averaged in the sense described 
above. Equation (25) is valid on time scales large compared to the period of 
quasioscillations of tp(t). This implies the following effective equation of state for 
the matter described by the scalar field (p 

P= q —^P- (26) 

q + 2 

In particular, for the most interesting cases q — 2 (quadratic potential) and q = 4 
(quadric potential) we obtain, respectively, p = (dust) and p = p/3 (radiation). 

Using the definition (15) we immediately obtain the approximate equation for 
the evolution of the value (po(t) 

00 « ^Hw (27) 



3. Particle production: Born approximation 

In this section we will study the process of particle production by the rapidly 
oscillating inflaton field. We take the interaction Lagrangian to be 



J int 



-f^ip - (a<p + hip 2 ) x 2 , (28) 



where ip and ip describe spinor particles and \ describes scalar particles with cor- 
responding masses mj, and m x , f and h are dimensionless coupling constants, and 
a is a coupling constant of dimension of mass. We will treat the scalar field p> as a 
classical external field. In this section we will employ ordinary perturbation theory 
in coupling constants, that is, we will work in the Born approximation. To calcu- 
late the rate of particle production we first develop the quasiperiodic evolution of 
the scalar field <p into harmonics 

00 

cp(t) = y^ (p n cos(nujt) , (29) 

71=1 

00 

<p 2 {t) = ^+^2( n cos{2nujt), (30) 

n=l 

where the value ip 2 m p>q/2 which is slowly varying with time is (p averaged over 
the rapid oscillations of the scalar field tp. <p n , ( n are amplitudes which are slowly 
varying with time, and 

uj = — = c^- (31) 

T <po 

is the leading frequency, also slowly varying with time, related to the period T of 

the oscillations of the field (p. The constant c in (31) is of order 1, and is given by 

-l 

in the case of a scalar field potential given by (2). To derive the estimate (31) we 
made use of the relation (14). 
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Assuming the condition (18), it follows that the oscillation period of tp is small 
compared to the Hubble time, 

H < uj (33) 

This last condition allows us to disregard the effect of cosmological expansion in 
evaluating particle production rates. We will also assume that particle masses and 
the coupling constant h are sufficiently small, so that 

m 2 , m 2 , + h^ < uj 2 . (34) 

The rates of particle production, that is, the total number of pairs produced 
per unit volume and unit time are then given in first order perturbation theory in 
coupling constants /, o and h by the following relations 



c2 x 



f 2 

H> = TO^ 5Z ^nuf , (35) 



n=l 



1 / CXJ 00 \ 

\ n—1 n—1 / 



Remembering the expansions (29) and (30) we can put the sums in the last ex- 
pressions into the following form 

oo 

J2 (fnnuj) 2 = 2 < 2 > , (37) 

n=l 



]>>2=2<y2 2 >, (3* 



71=1 



DC 



ECn = 2<(^ 2 -^ 2 ) >, (39) 

71=1 

where the brackets < . . . > denote the average over quasioscillations of the scalar 
field (fit). These average values in the expressions (37)- (39) can be estimated in 
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exactly the same manner in which the equation (24) is derived. We obtain, finally, 
the following relations for the particle production rates 



f q 

W M> = IZ^TTT.Pi (40) 



W 4vr q + 2 ' ' 



W 2 X = g^( c <r0"Vo +c/i/iVo) » ( 41 ) 



where 



■/ 

-1 


9 


dx 


if dx 


Vi 


-xi 


1 J y/l - xi ' 

-l 


f(x 


2 _ 


4)* 

-xi 


l 
dx j f dx 


1 

1 


1 J y/l-x9 

-1 



(42) 



Ch= / V K-Z^ -/ ^==, (43) 



are constants of order unity. 

Higher order contributions to the particle production rates (40)-(41) will be 
negligible if 

/ \ 2 



f, h, — « — . (44) 

Hence, our assumption (34) for the smallness of the coupling constant h is justified 
since (p 2 ~ y?o/2. We will see below what the conditions (44) imply for the 
concrete inflaton potentials. 

Equations (40)-(41) enable us to estimate the energy density production in 
the form of ip, ip and \ particles created. To do this, note that pairs of particles 
produced have average energy of the order u, if we reasonably assume that the 
sums on the right-hand sides of (35), (36) rapidly converge. Then the equation for 

12 



Pp, the particle energy density produced, can be written as follows 

p p (production) = (T^ + I\) p , (45) 

where 

I> = 0{l) f -u, (46) 

is the contribution corresponding to the production of fermionic particles, and 

T x = 0(1)^- + 0(1)^-, (47) 

X V ' 4:7100 v ' Attuj v ' 

is the contribution corresponding to the production of scalar particles \- Here and 
below we denote by the symbol 0(1) a constant of order unity. The values T^ 
and T x can be interpreted as the 'V-paxticle" decay rates in corresponding decay 
channels. 



4. Particle production: parametric resonance 

The expressions (45)-(47) for the particle energy density production rates have 
been obtained in first order Born approximation, and in this form they have been 
used in all the seminal papers on reheating. However, as it has been shown in [10], 
this approximation does not take appropriate account of the parametric resonance 
phenomena which take place during the oscillation period of the scalar field (p. 
Parametric resonance effects can be shown to be insignificant for fermionic fields 
[11]. The estimate (46) for fermion production made in the preceding section thus 
remains valid. However the question of the parametric resonance effect for the 
scalar field x has to be considered and this will be the aim of this chapter. We 
will show that the resonance occurs when the frequency oo^ of the quantum field 
mode is equal to the half- integer multiples of the inflaton frequency, uj\ rs (t| uj) . 
This results in exponentially enhanced production of particles in narrow resonance 
bands, at rates which are computed below. 

13 



The evolution of a particular mode Xk of the quantum scalar field \ i n the 
presence of the quickly oscillating classical scalar field tp with the interaction (28) 
is described by the following equation 

Xk + ZHxk + {k 2 + m\ + 2a V + 2h V 2 ) X k = , (48) 

where k = k/a is the physical wavenumber of the mode under consideration, k is 
the comoving number. 

Performing the transformation 

* = % m 

one obtains from the equation (48) the following equation for the function Y&: 

Y k +(ul(t)+g(ut))Y k = 0, (50) 



where 



.2/ + \ _ v2 , m 2 9 u2 3 - ~-^T 



and 



u&t) = k 2 + m 2 x - -H 2 --H + 2htp 2 , (51) 



g(wt) = 2a V + 2h Up 2 - <p 2 \ . (52) 



The function g(ut) is to a good approximation a 27r/w-periodic function of time t 
with two leading frequencies u and 2u>, as can be seen from (52) and the expansions 
(29), (30). We will regard this function as a small perturbation in the equation 
(50), and to be able to do this we require the condition 

2a(po + htfiQ ^C uj\ . (53) 

which we will assume to be valid. 
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The time dependence of the coefficients in the equations (48), (50) will lead to 
X-particle production which can be described by the standard Bogolyubov trans- 
formation technique. The problem consists in finding the solution to the equation 
(50) with initial conditions corresponding to the initial vacuum state and then 
evaluating the Bogolyubov transformation coefficients. We expect that for the 
equation (50) parametric resonance effect will be dominant hence we are going to 
study this effect in detail. 

From the general theory of parametric resonance (see Appendix A for details ) 
one knows that parametric resonance occurs for the equation (50) for certain values 
of the frequencies uu^. Namely, the resonance in the lowest frequency resonance 
band occurs for those values of oo^ for which 



u> 



2 tn \ 2 

k 



- l-u) = A n < \g n \, (54) 



2 J 



where n is an integer, and g n is the amplitude of the ra-th Fourier harmonic of 
the function g(ut). Of course, for such a conclusion to be valid in our case when 
the frequencies u^ and u depend on time we must be sure that these frequencies 
change with time slowly as compared to their own values, that is 



Uk 



LO k 



<^ ujk, w- (55) 



But this condition immediately follows from the estimates (20), (33), and from 
(53) if we also take the relation (54) into account. Note that condition (53) for uo k 
in the resonance band coincides with the conditions (44) on the validity of Born 
approximation. 

In this paper we will consider only resonance in the lowest resonance frequency 
band for each Fourier harmonic of the function g as given by equation (54). The 
effects of higher resonance bands, as is well known, are of higher order in the 



1) In Appendix A we explicitly introduced a small dimcnsionlcss parameter e. In the main 
text we have absorbed e into the definition of g. 
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amplitudes \g n \ which we assumed to be small. Hence we expect these effects to 
be less efficient. 

From the conditions (20), (33), (34) and (53) it also follows that part of the 
expression (51) satisfies 






u 2 k -k 2 = m 2 Y --H 2 --H + 2h<p 2 < uj 2 . (56) 



It then follows from the equation (54) that 

k 2 > ml - -H 2 --H + 21UP 2 (57) 

for the resonance values of k . This means that the redshift of the frequency u k 
due to the expansion of the universe is, to a good approximation, 

u k (t) oc k oc a _1 (t). (58) 

Due to the condition (33) we can neglect this cosmological redshift on the 
timescale of the scalar field tp oscllations. The mean occupation numbers for the 
generated x _ P a rticles are then given by the Bogolyubov coefficients /3& as JVfc = 
\Pk\ 2 an d the expression for them can be approximately written as follows (see 
Appendix B for derivation ) 

N k ~ sinh 2 [ I fi+dt) , (59) 



where 



1 



lluJ 



V+ = —\\9n\ 2 -^ 2 n (60) 



is the eigenvalue of the growing resonance mode, and A n is given by (54). 



1) We must use the expression (B32) for the number of particles produced, because in our case 
the value /i+ varies with time. 

1G 



The approximation (59) is valid only if we can treat the value ^u+ as adiabatic. 
For our purposes it is sufficient to consider values of [/,+ close to the resonance 
(when A n = 0) and the variation of //+ with time is mostly due to the expansion 
of the universe. Then the adiabaticity condition (see the equation (B31)) will be 
as follows 

fi+ > H. (61) 

As time goes on, the mode with the particular wavenumber k gets out of the 
resonance band, and the value 7V& ceases to grow. The reason for the frequency 
uj], (initially in resonance) to get out of the resonance band is that the redshift 
evolution (58) of the frequency Uk is different from the evolution of the leading fre- 
quency uj and hence the resonance frequency nuo/2 is continuously passing through 
the spectrum of frequencies u^ of the modes of the scalar field x- (I n the case 
of Ay? 4 potential, the value u does redshift at the same rate, and one has to take 
into account the back-reaction of the created particles to determine this frequency 
shift; see Section 5.) 

Denote by u Tes the resonance frequency nu/2 and by Au;& denote the difference 
uJk—uJros- Both frequencies u; res and u>k change with time. In a small time St the shift 
between these frequencies will be Sou = \dAuk/dt\ / \ UJk= Q5t. Hence a new region in 
phase space of volume u^ es 5uj/2ir 2 will be filled with particles. The average number 
of particles in every state in this phase space region will be given by N^ and each 
particle will have energy u; res . Then taking sum over all parametric resonance 
bands (with various numbers n in (54)) we can write down the equation for the 
energy density p x of the scalar particles produced in the following form 

p y (resonance production) = N — „ N TCS cjjL — - — , (62) 



2tt^ 

all resonance bands 



Au> fc =0 



where iV re s is the maximal value of JVfc, for the current resonance value of k, achieved 
after the corresponding mode has passed through the resonance band and has been 
amplified. 
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Under condition (34) particles produced will be ultrarelativistic and hence their 
contribution to energy density will decrease due to the cosmological expansion. 
Taking also this process into account we are able to write down the following com- 
plete equation for the evolution of the energy density p p of the particles produced 
(including the contribution from fermions) 

Pp = -ah Pp + (r^ + ri res) ) P , (63) 



where 



r (res) _ 1 V^ N r , 3 

X ~ 2n 2 p ^ JV res^ rC s 

all resonance bands 



dAook 



(It 



(64) 

Auj k =0 



plays the role analogous to that of T x in (45). The values of r^ given by (64) and 
T x given by (47) are to be compared in order to conclude about the significance 
of the parametric resonance for the process of scalar particle production. In what 
follows we consider two important cases of the inflaton potential V((p), namely, 
quartic and quadratic potentials. 



5. Case of potential V((p) = \(f 4 
In this case we have from (31) 

u> = c v\<po , (65) 

where the constant c is given by (32) with q — 4. The resonance frequency u rca 
is thus proportional to tpo and by virtue of (27) in the absence of back reaction of 
the particle production process on the scalar field (p we would have 

u Tea oc a~ . (66) 

As we also have the condition (58) this would mean that the resonance amplification 
takes place for all times in the same wavebands, which satisfy the condition (54). 
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However due to the energy loss due to particle production, the scalar field amplitude 
tpo, and hence the frequency u; res , will decrease more rapidly than given by (66). 
Wee see that in order to make correct estimates of particle production in the model 
with quartic inflaton potential one has to take into account back-reaction on the 
inflaton evolution. 

Let us then consider back-reaction of the particle production process on the 
evolution of the scalar field tp. Since particle production is a random quantum 
mechanical phenomenon, the change in the evolution of the scalar field tp will be 
random. We can picture the back-reaction as leading to abrupt changes in the 
state of the scalar field at those random moments of time at which particles are 
produced. On sufficiently large time scales, on which many particles are produced, 
this change will exhibit regular features, and it is these features which we are to 
determine. 

The evolution of the scalar field tp can be described by two parameters, namely, 
the scalar field energy density p and the frequency uj of its quasioscillations. In 
absence of back-reaction both these parameters are changing slowly with time, 
when compared to the time scale T = 2tt/uj of the quasioscillations. The value u is 
directly related to the energy density through the relation (31). We find it possible 
to describe the state of the scalar field in terms of these two parameters also if we 
take the back-reaction mentioned above into account. Indeed, it is always possible 
to speak of the scalar field energy density. As far as its quasioscillation frequency is 
concerned we note that its very existence is necessary for particle production, with 
rates given by (40)-(41), or with average numbers (59) in the case of parametric 
resonance. 

Also note that on the quasioscillation time scale there can be but few par- 
ticles produced (this fairly general statement follows from the standard analysis 
of transition amplitudes). Hence random moments of the abrupt changes in the 
scalar field state will be time-separated in average by more than one scalar field 
quasioscillation period. Between these moments of abrupt changes the scalar field 
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is evolving unperturbed. From this it follows that we can assume the validity of 
the relation (31) between the quasioscillation frequency uo and the energy density 
p. 

To determine in the most simple way the changes in the evolution of the scalar 
field energy density p we can apply the energy conservation law. The energy is 
transferred from the scalar field (p to the particles produced. Hence, given the 
energy production rate (63) of the latter we thereby know the energy loss of the 
scalar field (p. The equation for the inflaton energy density will then be 



^(res) 



p = -(4i7 + r^ + rnp, (67) 



^(res) . 



with T^ given by (64). From the equation (2) with q = 4 and from (15) it then 
follows 

00 = - (h + \ (i> + rjf"*)) <p . (68) 

Now taking (58) and (65) into account we will obtain the following expression for 
the time derivative factor in (62) 



dAuk 



(It 



Aw fc =0 4 



l(r> + r£ es) ) tu ms , (69) 



and from the equation (64) we then obtain 



J2 N res u;i s = 8n 2 p[l + ^-) . (70) 




all resonance bands 



It is very remarkable that in spite of all the details of the interaction between 
the fields (p and x the maximal occupation number mean values depend through 
(70) only on the resonance frequencies u; res in different resonance bands. This 
significant feature of the quartic inflaton potential makes the total analysis rather 
simple. 
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It remains to estimate the value 1^ . The expression for 7V& is given by (59) 
where the integration in the argument of the hyperbolic sine is taken over the time 
during which resonance is taking place for the mode with the wavenumber k. In 
terms of the value Au^ = oo k — u; res , the expression (60) for //+ can be rewritten 
as follows 

1 



H+ = 7, V \9n\ 2 - Au 2 k (Au k + 2u res ) 2 . (71) 

Then using (69) and taking into account the smallness of the value \g n \ as compared 
to uj 2 es (this follows from the relations (52) and (53)) we can estimate the integral 
in the exponent of (59) as follows 

\g„\/2LOrcs 

1 



lX + dt « / 7TT — vlS'ril 2 — A^|^2u; rcs )2 _ __ dAbJ k 

— \g„\/2oj T , 



^I'C'.T 



2 ^ es v— - --- ^ + r^,^, 



2 r m- l„ |2 

1 — x 2 dx 



\9n\ I _ A 2 J- n \9n 



r^ + r^u^l z^ + vt^L 

(72) 
In deriving (72) we considered slowly changing variables as constants. Then the 
expression (59) for N res becomes 



Nres * sinh 2 ( * l9 f I . (73) 



Substituting this expression for N res into (70) we obtain the following implicit 
expression for the value T^ 



u ™ smh —, 7— rr ~ 87T p 1 + -^ 

l 2(r. + rt res) U r 3 / v r («) 



|2 

■es 
all resonance bands \^l- L V'~'~''"X J ^res / V i X 

(74) 
Consider first the case when energy production of scalar particles dominates 
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over that of fermions so that 



r> < it es) . (75) 



Typically one of the terms, say, of the n-th resonance band dominates in the sum 
in the left-hand side of (74). In the case (75) 1^ is determined by this dominant 
term as follows 

I*-) K *£ log- (BpL) . (76) 



^res V ^res 



Two of the Fourier harmonics of the function g given by (52) will be of appre- 
ciable value, namely, 



\gi\ ~ <T(po, (with u res = —J , (77) 



and 



\g 2 \ ~ — ^, (with u; res = u) . (78) 



If the contribution from n — 1 is dominant in (74) (that is if (po < cr/h) then 
from (76), (77) and using the expression (65) for uo and the explicit expression for 
the inflaton energy density p = Ay?Q we get 



,(res) 8tt a 2 _j /512tt 2 \ a 



In this case the adiabaticity condition (61) reads 



9 Mpa 

tpl « -^-. (80) 



If the contribution from n = 2 is dominant in (74) then from (76), (78) we in 
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a similar way obtain 



(res) nh 2 -i/32tt 2 \ 

r x -1^372^0 log (^j, forv, > 



In this case the adiabaticity condition (61) implies 



ipo < M P -. (82) 



According to the note made at the end of the previous section, the values (79) 
and (81) for T^ are to be compared to the values of T x in the corresponding 
cases. The expression for T x is given by (47). Calculating the ratios of T x to T x 
in the two cases discussed just above we obtain 



r (res) 
L X 

Tx 


r^j 


32 7T 2 

c 2 A 


log" 1 | 


/512tt 2 \ 

UaJ' 


for (po < - , 
h 


r (res) 
L X 

Tx 


rs_, 


7T 2 

c 2 X 


log" 1 | 


/32tt 2 \ 

U 4 aJ' 


for ip > - . 
h 
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We see that these ratios look similar in two cases considered, and for typical 
values of the inflaton coupling constant [1] 

A ~ 1(T 12 (85) 

they are extremely large. Hence, for the inflaton field with the potential A<^ 4 with 
small A, in the domain of applicability of perturbation theory (given by (44)), 
parametric resonance effects strongly dominate over the usual Born approximation 
estimates, so that the latter represent a serious underestimation of the particle 
production rate. 
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Let us proceed further and estimate the possible reheating temperature for the 
model considered. The equation for the energy density of the particles produced 
is (63). We should take into account also the equation (67), modify the equation 
(12) as follows 



a _ 8tt 
3Af2 



^ = ^2 (p + Pp)> (86) 



and then solve the system of equations (63), (67) and (86) with initial condition 
p p = at the moment of the end of inflation. Then the maximal value achieved 
in the course of evolution by p p would determine the reheating temperature. The 
system of equations mentioned can be solved exactly if we use the condition (75) to 
neglect T^ and the expressions (79) or (81) for 1^ . In the case (79) the solution 
will be 

'^(/-i^- 1 )) 4 ^ (87) 

and in the case (81) the solution is 



»i 1+ m leT - 1) )~ i ''~ ir - 



where To, and Hq are the initial values of T^ r and H respectively, po and ao are 
the initial values of the total energy density and the scale factor, and r = log (a/ao). 
In both cases (87) and (88) the solution for p p is given by 



P P = Poe- 4r -p- (89) 



Until the condition p p ^ p holds, the solutions (87) and (88) can be approxi- 



1) We take the reheating temperature to be the temperature which would correspond to the 
energy density of the particles produced. This assumes that the particles produced will be 
thermalized soon after their creation. In this paper we do not touch on the theory of the 
thcrmalization. 
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mately written as 




(90) 



and 

o^H^-^r^' <9i) 

respectively. 

The solutions obtained determine the maximal value achieved by p p . Using 
the expressions (89), (90) and (91) it is easy to see that this value is achieved as 
soon as the condition 

r^ res) ;> AH (92) 

becomes valid. After this particle creation proceeds very effectively and the inflaton 
field loses all its energy in less then one Hubble time. We shall therefore use the 
condition (92) to estimate the reheating temperature in various cases. First, assume 
that reheating takes place in the region (p$ > o jh. Then using the expression (81) 
for T^ r we obtain from (92) 

A ) l ° gl A ' (93) 

where we also used that A is small to simplify the logarithmic factor. Now, the 
condition y?o > cr/h assumed above will read 

a & KT 1 M P h (^\ log" 1 -j- . (94) 

Hence if the condition (94) is valid, the value of the inflaton field amplitude <po 
at the moment of reheating is given by the right-hand side of (93). Note that the 
adiabaticity condition (82) is valid for the values of (/?o given by (93) since the ratio 
h/X has to be small according to the conditions (44). 
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If the condition (94) is not valid then reheating takes place in the region </?o < 
a/h and we must use the expression (79) for 1^ . In this case from (92) we obtain 
the estimate 

^ ~ A 2/3 lQ S 7 J ' (95) 

and the condition (p$ < a/h would imply, as it should, the opposite of (94), namely 

a > 10" 1 M P h(^\ log" 1 \- (96) 



In this case the last inequality together with (95) leads to the estimate 

2 



V? £ lO^Mp^ log" 1 ! 



(97) 



for the boundary value of </?o at which effective reheating starts. 

Due to the conditions (44) the adiabaticity condition (80) can again be seen 
to be valid for the values (95) of ipo. 

We see that the expression in the right-hand side of (93) gives the minimum 

possible value for the inflaton field amplitude (po at the moment of reheating. We 

stress again that for the values of </?o obeying (93) the condition AH ^ I* 

holds and reheating is extremely effective. To estimate the numerical value of the 

threshold (93) note, that according to the conditions (44) (which we remember 

to be necessary in order that perturbation theory is valid) we have the condition 

/i< A. So taking for example 

h i , , 

- ~ KT 1 (98) 

A 

and taking into account (85) we will get from (93) 

ifo & l(T 4 Mp. (99) 

The reheating temperature T r h is then estimated from the condition that the value 
of p p at the moment of reheating is of the order of the value of p at the moment 

2G 



when the condition (92) is first achieved, and from the relation 

T rh ~ pl ,A . (100) 

Finally we obtain 

T rh ~ A 1/4 y?o ~ 10" 1 M P A 1/4 f^\ log" 1 i (101) 

In case of the values (85), (98) for the couplings we get 

T rh ~ 10~ 7 M P ~ 10 12 GeV. (102) 

Now let us see what our condition (75) means. In the region </?o > &/h using 
the expressions (46) and (81) we obtain the condition for the coupling / (again 
simplifying logarithmic factor by using smallness of the coupling constant A) 

f 2 < O(10) (j^J log" 1 1 , (103) 

which does not depend on the value of </?o- For the values (98) of the ratio h/X the 
condition (103) is not very restrictive. 

In the region </?o < &/h using expression (79) we obtain from (75) the estimate 

which depends on the value of y?o but certainly is satisfied for (po sufficiently small. 
For the condition (104) to be valid for our estimate (95) of the inflaton field am- 
plitude at the moment of reheating the coupling / is to be 

f 2 <°^(jk) log " I/3 r (105) 

Thus, to conclude this part, if the value a satisfies (94) and the coupling / 
satisfies (103) then the value of the inflaton field amplitude y?o at the moment of 
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reheating is given by the right-hand side of (93). If the value a satisfies (96) and 
the coupling / satisfies (105) then the value of the inflaton field amplitude <po at 
the moment of reheating is given by the right-hand side of (95) and satisfies (97). 
In both cases the reheating temperature is not less then given by (101), with the 
numerical estimates (102) for the coupling values (85) and (98). 

If neither of the two mutual conditions just mentioned holds then at the mo- 
ment of reheating the condition (75) is not valid, and the reheating temperature 
is determined by fermionic particle production. In this latter case at the moment 
of reheating we have 

ri rcs) < i>, (106) 

and in the condition (92) for reheating we should replace the value T-^ by T^. 
Then using expression (46) for T^ we obtain for the value of the inflaton field 
amplitude </?o at the moment of effective reheating the well-known expression 

V? £ 10- 2 / 2 M P , (107) 

and for the reheating temperature we get 

T rh ~ 10~ 2 A 1/4 f M P . (108) 



In the case (106) from the equation (74) we obtain the following expression for 



r (res) 



\'2 
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It is useful to compare this value of T^ to the value of T x given by (47). To 
do this note that mainly two resonance bands contribute to the sum (109), namely, 
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— sinh 




r„ 




2-71 


^n 
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K\g n \ 2 


2X ^ ui Tes 



those given by (77) and (78). Therefore we can write 

r (res) = r (res) (cr)+r (rcs) (/ ^ (nQ) 

where T\ depends only on the coupling a, and T^ only on the coupling h. 
These two terms in the sum (110) correspond to the two terms in the expression 
(47) for T x . We then rewrite the expression for T x in the similar way 

r x = ri((7) + r 2 (/o. (in) 

Now using the expressions (46), (47), (65), (77) and (78) we can easily estimate 
the ratios of each of the two terms in (110) to its corresponding term in (111) as 



(112) 

where 

(113) 

We can see that parametric resonance effects for the scalar field \ would be 
significant in the case (106) if the values of x n given by (113) were larger then unity. 
This last condition however can be shown to be incompatible with the assumption 
that fermionic particle production dominates at the moment of efficient reheating. 
Hence at this moment the value T x given by (47) is larger than the resonance value 
1^ given under the conditions considered by (109). 

All the results obtained in this section for the inflaton field amplitude during 
effective reheating can be written in a single equation 

** s Mp max { 1(r2/2 ' 10 " (ff log " v dhf' 3 

' :n4) 

from which also all the necessary inequalities between coupling constants can be 
easily derived. We stress once again that for the values of (po given by (114) the 
condition AH < T^ + T x is valid so the reheating process is very effective. 
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The expression for the temperature of reheating is then 



M P AV 4max { 10 - 2/ M0-(^) 2 lo g -i,(^ 



2/3 , -1/3 ( 115 ) 

log l ' 6 - 



Mp\J ° A 



It is very instructive to compare our results (114), (115) to those which would 
be obtained by making use of Born approximation (described in Section 3) for the 
scalar particles \ production rates. If instead of Y x we used the value T x given 
by (47) then we would arrive at the following reheating conditions 

V? £ M P max <j 1(T 2 / 2 , 10" 2 A (J\ , A (j^) } , (116) 



and 



rh 



Mp^ m ax{l0-V 2 , lO-A^) 2 , > (j^f} - (117) 



One can see that for the scalar particle couplings the values (116), (117) are smaller 
by an enormous factor [A log(l/A)] than the corresponding correct estimates 
(114) and (115). 



1 ™2,^2 



6. Case of potential V(cp) = ^m (p 

In this section the analysis will procede along the same lines as in the previous 
one. In the case of a quadratic inflaton potential we get from (31), (32) 



118) 



so the inflaton field oscillates with constant frequency. The resonance frequency 
o; r es = nuo/2 will thus be constant in time. For the time derivative factor in (62) 
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we then obtain, taking into account (58 

dAuk 



dt 



= H UJ lea , (119) 

Au> fe =0 



and from (64) we will have 

rt S) = ^- ^ N res u« es , (120) 



,(res) _ H _ S^ AT i ,4 

' P 
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where all the relevant variables have been explained in Section 4. The expressions 
(59), (71) remain valid in the case considered but the estimate analogous to (72) 
will be modified as follows 

\g n \/2uj vca 

/ fx+dt « / J\g n \ 2 - Au;|(2u; rGS ) 2 — dAu k « 

7 J *U} Tes v tl a; re s 



|Sn|/2W: 



res 



:i2r 



i 

l^l 2 /" A 2^ ^ l^l 2 

1 — x A dx = — 



4tfu; r 3 es 7 • ' 8iJa; r 3 es ' 

-l 



so that the occupation numbers iV res will be given by 



.V,,, ~ sm h 2 f?#£-l . (122) 



8 Hul 



res 



^(res) 



and the equation (120) for 1^ will read 

r (res) ^ _B_ ^ 4 2 /tt \ 9n \ 2 \ 



(123) 
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We must compare this value to the value of T x obtained in the Born approxima- 
tion. Writing T^ and T x in the form (110) and (111) and using the expressions 
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(47), (77) and (78), we will arrive at an estimate of the same form as (112) but 
with different x n : 

I |2 
71 \9n\ /io/i\ 

8 #^ es 
Again we see that parametric resonance effects for the scalar field x are significant 
if the values x n are large. 

As usual, one of the terms in the sum (123) will dominate. In the case </?o < cr/h 
the harmonic with n — 1 dominates, whereas for tpQ > a/h the harmonic corre- 
sponding to n = 2 is more important. The relevant expressions for the amplitudes 
in these two cases are given by (77) and (78). 

In the case tpo < a/h we get 

(res) „ Hm 2 h 2( ™^Po \ (m) 

ix " 167rVf \Hm*)' [ ' 

and the condition of effective particle production (92) in this case will read 

y?0\ 2 . Hm Hm 

— ) ^ o lo §^- ( 126 ) 

As soon as particle creation begins the Hubble parameter depends on the total 
energy density of the inflaton field and of the particles created. So the value of 
the Hubble parameter definitely exceeds the expression given by (12) with only 
inflaton energy density taken into account. This enables us to obtain the lowest 
possible bound on the region of the effective particle production. If we substitute 
in (126) the lowest possible value for the Hubble parameter given by (12) we will 
get the result, up to numerical factor of order one, 

— k, -2-rp- log -2-rr- • 127 

m o A Alp a A Mp 

This estimate is valid only if the coupling constant a is in the range given by 

< I — 1 < t^- (128) 



M P J \mJ M P 

The first inequality in (128) is the condition for perturbation theory to be valid, it 
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stems from (44). The second inequality in (128) is the criterion of large occupation 
numbers of the particles produced, that is, of effectiveness of the resonance particle 
production (in other terms, it is the condition for the value x\ given by (124) to 
be larger then unity). The validity of the adiabaticity condition (61) in this case 
follows from (127) and from the first inequality in (128). 

In any case there is also another region in which the condition (92) holds so 
that particle production can be effective. This region is given by 



y?0 < o l M v 
m mr 



(129) 



and it can be obtained by using the first term in the expression (47) for the value 
T x since parametric resonance in this region is not effective. 

In the case opposite to the one just considered, namely, when ipo > a/h we 
obtain 

and the condition (92) will give the estimate 

y?o\ 4 32H H 

— ) ~ ~n~ lo si2— > 131 

The lowest possible bound on the effective particle production region in this case 
will be obtained if we substitute the lowest possible value for H, given by (12), 
into the right-hand side. We will get 

y?o\ 3 . 65m m 

m) * FMp l0g FMp- (132) 

The theory which led to the result (132) will be self-consistent if the coupling h is 
in the range 

m \ , 2 . . m 



w i < h < T^ ■ ( 133 ) 

Again, the first inequality in (133) is the requirement that perturbation theory be 
valid and the second one is the requirement of large occupation numbers of the 
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particles produced. For the values of <po (132) the adiabaticity condition (61) can 
be shown to hold. 

We emphasize that the condition (92) in the case considered determines not 
the possible reheating temperature but rather the "freeze-out" boundary (127) or 
(132) for the inflaton field tp: below the values of <po given by these expressions the 
particle creation process is ineffective and the inflaton field energy density decrease 
is dominated by the cosmological redshift rather then by particle production. The 
energy density of ultrarelativistic particles created decreases more rapidly than 
that of the inflaton field and the universe can well become again dominated by the 
inflaton field (p until the condition (129) is satisfied. If the couplings a and h are 
extremely small so that 

then the effective reheating condition (92) in the model with quadratic inflaton 
potential takes place only for the values of tpo given by (129). It is easy to see that 
in all the cases considered above the minimal possible "freeze-out" energy density 
Pf reeze of the inflaton field is 

p freeze ~ m 4 . (135) 

If during the particle creation process the energy density of the universe be- 
comes dominated by the particles produced, then the estimates for the actual values 
of the inflaton field energy density at which it freezes out will be higher then those 
given by (127) and (132). Such a condition is likely to take place in the range of 
coupling constants given by (128) and (133) since, as discussed in [13], at the onset 
of the oscillation period of the inflaton, reheating process may be very efficient. In 
this case one must use the equations (126) and (131) in which it should be taken 
into account that the Hubble parameter H is dominated by the contribution from 
the particles created and, hence, is much higher then given by Eq. (18). We note 
that in the range of couplings (128) and (133), at the moment of the onset of the 
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inflaton oscillations, perturbation theory developed in this paper breaks down. It 
becomes valid again when the amplitude y?o decreases sufficiently and the condi- 
tions (44) start to hold. The value of the Hubble parameter in the equations (126) 
and (131), however, cannot be higher then its value H ~ m at the end of inflation. 
This allows one to estimate also the maximal possible boundaries for the inflaton 
field freezing out which will be [13] 

(^) 2 ~ ("flog™, (136) 

and 

(^Y ~ h- 2 \ogh-\ (137) 

\m / 

These estimates are just on the boundary of the applicability of perturbation the- 
ory. 

If interactions of the inflaton <p with fermions take place, we obtain, using the 
equation (46), the following estimate 

V?o £ 10~ 2 f 2 M P (138) 

for the reheating condition, which coinsides with the analogous estimate (107) in 
the case of the potential V((p) = Ay? 4 which we considered in the previous section. 
It can be shown, using the expression (31) for the value u and the expression (46) 
for the value of T^, that the estimate (138) does not depend on the power q of the 
inflaton potential V(ip). 
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7. Discussion 

In this paper we considered the problem of scalar particle production after in- 
flation by a quickly oscillating inflaton field. We were using the framework of the 
chaotic inflation scenario with quartic and quadratic inflaton potentials, and we 
considered linear and quadratic coupling of the inflaton field tp to a quantum scalar 
field x ( see formula (28) for the interaction Lagrangian). The particle production 
effect has been studied using the standard technique of Bogolyubov transforma- 
tions. Specific attention has been paid to parametric resonance phenomena which 
take place in the presence of a quickly oscillating inflaton field. 

We have found that in the region of applicability of perturbation theory (when 
inequalities (44) are valid) the effects of parametric resonance are crucial, and 
estimates based on first order Born approximation are often not correct. 

In the case of the quartic inflaton potential V((p) = Ay? 4 the reheating process 
is very efficient for any type of coupling between the inflaton field and the scalar 
field x even for small values of coupling constants. The reheating temperature in 
this case is [A log (1/A)]~ times larger than the corresponding estimates based on 
first order Born approximation (see the expressions (115), (117)). 

In the case of a quadratic inflaton potential the situation is more complicated. 
The reheating process depends crucially on the type of coupling between the infla- 
ton and the scalar field x an d on the magnitudes of the coupling constants a and 
h. The theory predicts not ony the possible reheating temperature but also the 
effective energy density of the inflaton field "freezing out" below which the infla- 
ton field energy density decrease is dominated by the cosmological redshift rather 
then by particle production. The inflaton, so to say, decouples from the matter 
it creates. In the case of the quadratic inflaton potential we obtained the low- 
est (equations (127) and (132)) and highest (equations (136) and (137)) possible 
boundaries of the effective energy density of the inflaton field when it freezes out. 
With an interaction linear in the inflaton field, besides these possible freeze-out 
boundaries which exist for couplings in the range (128), there is also a reheating 
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boundary (129). For an interaction quadratic in the inflaton field the lowest possi- 
ble freeze-out boundary (132) exists in the range (133) of the values of coupling h. 
If the couplings a and h are as small as to satisfy (134) then the effective reheating 
condition V ^ AH is achieved only at the later stages after the condition (129) 
starts to hold. 

The situation is more complex for more complicated potentials. For example, 
one may wish to consider the inflaton potential of shape 

V(<p) = \{ V 2 -rf)\ (139) 

with spontaneous symmetry breaking at a scale 77. In the chaotic inflation scenario 
the inflaton field rolls from large values of \<p\ to the minimum of its potential. 
After inflation it starts quasiperiodic motion. For values of the amplitude y?o ^ V 
the behaviour of the inflaton is like for a Ay? 4 potential, and reheating proceeds 
as described in Section 5. But as soon as the magnitude of y?o becomes close to 
77 the dynamics changes and the inflaton field starts oscillating around one of the 
minima of its potential at \tp\ = 77. The inflaton potential in the vicinity of its 
minimum can be approximated as quadratic, hence the reheating (freeze-out) will 
be described by the theory of Section 6. 

In parallel to our own work, Kofman, Linde and Starobinsky have also been 
working on these issues for several years. Their results were recently summarized 
in Ref. 13. While our paper deals only with the case when the conditions (44) 
hold and perturbation theory in the inflaton couplings to the other fields is valid, 
both the cases (44) and opposite to (44) are considered in [13]. The resonance 
frequency bandwidth can be shown to be relatively narrow in the case (44), and 
broad in the case opposite to (44), so these two cases are called in [13] "narrow" 
and "broad" resonance cases respectively. In the case of broad resonance, higher 
resonance bands dominate the contribution to the particle production. After the 
inflaton field amplitude decreases sufficiently one enters into the regime of narrow 
resonance which can be accurately described by perturbation theory. The results 
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announced in [13] for the case of narrow resonance in general agree with those 
obtained in our paper in Chapters 5 and 6. There are some results for the case of 
broad resonance which we would like to compare with our results for the narrow 
resonance case. In doing this we shall use the notation of the present paper which 
somewhat differs from that of [13]. 

In the case of quartic inflaton potential the estimate of [13] Nf*(t) oc exp ( ^vA^omI 
for the case of coupling h = 6A (self-excitation of the inflaton, considered in [13]) 
can be compared with our estimate Nk(t) oc exp(2//+£) ~ exp ( -^yXipot) for the 
case of smaller coupling (/i<C A). We note that these estimates will be numerically 
of the same order of magnitude if we choose our smaller coupling to be h ~ 2A/5. 
For the values of h ~ 0.22A our estimate (93) for the beginning of reheating will 
be of the same order as the estimate of [13] y?o ^ 0.005Mp log" (1/A) for h = 6A. 
These estimates show that the effect of other particle production for the values of 
coupling about h <; 0.3A may be larger than the effect of inflaton self-excitation. 
This value of h is on the border of the applicability of perturbation theory. We 
cannot say anything about the reheating with couplings h ^ A as our perturbation 
theory breaks down in this case. 

If the inflaton potential is quadratic then, as previously discussed by Kofman, 
Linde and Starobinsky (see Ref. 13), in the case of absence of inflaton coupling 
to fermions and without linear (in inflaton field) coupling to scalar particles, the 
inflaton field will eventually decouple from the rest of the matter, and the residual 
inflaton oscillations may provide the (cold) dark matter of the universe. Decoupling 
of the inflaton field occurs somewhere between the values given by the right-hand 
sides of (132) (lowest possible freeze-out boundary) and (137) (highest possible 
freeze-out boundary). 
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Appendix A 

In this appendix we present necessary expressions which describe the effect of 
parametric resonance in the lowest resonance band. 

Consider the equation for the function Y(t) of the following form 

Y + (cut + eg(wt)) Y = 0, (Al) 

where uq and u are constant parameters, g(x) is a 27r-periodic function, and e is a 
small number. Our task is to find the approximate solution to the equation (Al) in 
the first instability band of the frequencies c^o- To do this we will use Bogolyubov 
method of averaging which implies perturbation theory of certain type in the 
small parameter e. 

The function g(x) being periodic, it can be developed in Fourier series as follows 

00 

g(x) = Y, 9ne inx , (A2) 

71= — OO 

with the amplitudes g n satisfying 

g*n = 9-n- (A3) 

With no loss of generality (for a small parameter e) we can put go = (redefining 
c^o if necessary). For the following purposes it will be convenient to introduce the 
phases a n as follows 

g n = \g n \e ia ". (A4) 



^ = (P u ) +e A, (A5) 



Let us introduce the value A by the following relation 

2 

where p/q is a rational non-contractible number (p and q are integers). The equa 
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tion (Al) then becomes 

Y+(^uA Y = -e(g(ut) + A)Y. (A6) 

The Bogolyubov method [14] consists in looking for a solution to the equation 
(A6) in the form of the following expansion in powers of e 

Y = acos^ + JV« (s) (a, 9, -t) , (A7) 



where 



.=1 v 9 



1/; = lujt + 9, (A8) 

Q 



the u^ s ' are functions periodic in their second and third arguments, a is the am- 
plitude and 9 is the phase of the solution (A7). The values of a and 9 are not 
constant in the Bogolyubov approach. Rather, they are functions of time t. 

We will be concerned explicitly only with the first order approximation in the 
small parameter e. For the function u^- 1 ' we have 



t (1) f o, 9, -t) = ^2 u n(a, 



«W o, 0, -tj = >^ Mn(a^)e ! V, (A9) 

with coefficients u n (a, 9) periodic in 9. 

Clearly, in the case of e = the solution to (A6) will be just the first term in 

the right-hand side of (A7), with arbitrary constant amplitude a and phase 9, so in 

this case we have a = 0, 9 = 0. For e^ according to the Bogolyubov method we 

regard the values d and 9 as functions of a and 9 and we expand these functions 

in powers of e 

a = eA(a, 6) + O (e 2 ) , 

(A10) 
9 = eB(a, 9) + O (e 2 ) . 

Putting the anzatz (A7) into the equation (A6) and using (A10), after collecting 
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terms of first order in e we will get the following equation for the function u^- 1 ' 
u^+I-lj) m (1) = -(^ + A)acos^ + 2-a;i4siii^ + 2-a;Socos^. (All) 



Using the expansions (A2) and (A9) we can develop both sides of the equation 
(All) in harmonics exp(i^ut) . On the left-hand side of (All) we will have, to the 
zeroth order in e, 

n——oo ^™ ' \y/ 

and we see that terms with n = ±p in this expansion are equal to zero. The 
corresponding terms in the development of the right-hand side of (All) should 
also vanish. These conditions allow us to determine the functions A(a, 9) and 
B(a, 6). If we denote by s the number 



2p 
s = — , 
Q 

then the expressions for these functions will look as follows 



(A13) 



A(a, 6) = ^ sin(2# - a s ) , 

sou 



B(a,6) = —(A + \g $ \cos(29-a s )) 

SUJ 



(A14) 



with the phases a s defined in (A4) for s integer. The expressions (A14) are valid 
for arbitrary s, in the case of s non-integer we should simply put g s to zero, and 
the (unspecified) value of a s is of no importance. 

If we make the change of variables 

fa as \ 
x = a cos I U —I , 

(A15) 

y = a sin [9 - y^ 
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then the system (A10) in terms of new variables x and y will look as follows 

x = —(\g s \-A)y + 0(e 2 ) , 

SCO 

(A16) 
y = ±(\g s \ + A)x + 0(e 2 ) . 



SijJ 



To first order in e the system (A16) is linear and its eigenvalues are 



pt ± = ± — J\g,\*-tf. (A17) 

suj v 

From this expression one can see that fi± can be real only if g s ^ 0, that is if s is 
integer. Hence unstable growth of the solution Y to the equation (A6) can take 
place only for the values of A (A is defined in (A5)) 

A < \g a \, (A18) 

and the growth (decay) rate is exponential, Y ~ exp(/i±t), with /i± given by 
(A17). 

For every integer s for which g s is non-zero the expression (A17) gives the 
growth rate of the solution Y in the first instability band determined by (A5) 
and (A18). There are also infinite number of instability bands whose width and 
growth rates are higher order in e (to all orders of e). Relevant expressions can be 
obtained to any desirable order in e by keeping track of the higher order terms in 
the expansions (A7) and (A10). For the purpose of this paper these higher order 
expressions are not required. 

Appendix B 

The aim of this Appendix is to provide necessary detailes of the resonance par- 
ticle production effect studied in Section 3 of this paper, in particular, to derive the 
formula (60). We first consider the general problem of a harmonic oscillator with 
time-dependent frequency, and after that the problem of resonant scalar particle 
production. 
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a) Harmonic oscillator with time-dependent frequency. 

Consider an oscillator with coordinate Q, conjugate momentum P, whose fre- 
quency Q depends on time so that the Hamiltonian is 

n = \ (p 2 + n 2 (t)Q 2 ) . (bi) 

For a quantum oscillator the values Q and P are Hermitian operators with 
standard commutation relations (we put % — 1) 

[Q,P] = i. (B2) 

In the Schrodinger representation, these operators are time-independent. 

Our aim will be to go to a convenient time-dependent frame in Hilbert space, 
in which the Hamiltonian (Bl) is diagonal at every moment of time. To do this 
define time-dependent operators a(t) and a*(t) by 

i J Qdt 

{QQ + iP) , 



2Q 

(B3) 



-ijQdt 

a) = =^- (VLQ - iP) . 



2Q 

These operators are mutually Hermitian conjugate and have the standard commu- 
tation relations of creation-annihilation operators 

[a, a f ] = 1 . (B4) 

In terms of these operators the Hamiltonian is expressed as follows 

H = tt(- + a^a] . (B5) 

The orthonormal frame in Hilbert space which diagonalizes the Hamiltonian 
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H is given by time-dependent states 



\n t > = K )4 \0t > , (B6) 

vn! 



constructed from the time- dependent vacuum state |0j > which is annihilated by 
the operator a(t). 

The operators a(t) and a' (£) obey the following equations of motion 



n 

(B7) 



_ e 2ifSldt f 

2fi 



a 1 = — — ■ e "" ■> """"" a . 



,t — „-2i ffidt 



so these operators vary with time only when Q varies . 

The solution to the equations (B7) can be written in terms of constant creation 
and annihilation operators aj and an as follows 

ait) = a(t) an + 8* (t) al , 

o ^ 

a\t) = f3(t) an + a* (t) a\ , 

with the complex functions a[t) and 8(t) obeying the system of equations similar 
to (B7) 

n 



2Q ' 



P = H e -2iJVdt a ^ 

1 2Q 



(B9) 



The system (B8) represents what is called Bogolyubov transformation between two 
pairs of creation - annihilation operators. 



1) This property is achieved by a convenient choice of the time-dependent phase factors in the 
definition (B3) of the operators a and a' . 
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If the oscillator initially (at t — 0) is in the vacuum state then its state |Oo > 
is annihilated by the operator ao and the initial conditions for the functions a(t) 
and f3{t) are 

| a (0)| = 1 , /3(0) = . (BIO) 

At the moment t it will not be in the vacuum state \0t >, annihilated by the 
operator a(t). Rather, there is the following relation between the states considered 

|0 > = ; l exp (^$- (aUt)) 2 ] 10* > , (Bll) 

which follows from the equations (B8). At the moment t the average number of 
the oscillator excitation level (number of quanta) is 

N(t) =<0 |a t (t)a(t)|0 >= \P(t)\ 2 ■ (B12) 

In the Heisenberg representation, the equation of motion for the operator Q(t) 

is 

Q + n 2 Q = 0. (B13) 

The solution to this equation can be expressed in terms of the operators ao an d % 
so that 

Q(t) = Q^-\t)a + Q {+ \t)al. (B14) 

Comparing this expression with the equations (B3) and (B8) we can express the 
coefficient (5{t) in terms of Q^'(t) at the moment of time at which the first time 
derivative of the frequency Q vanishes as follows 

13 = ^=e- { $ mt foQ H - iQ { - ] ) . (B15) 



2n 

b) Resonance particle production. 

4G 



Now let us turn to the issue of scalar particle creation in the external spatially 
homogeneous periodic field eg (out) whose properties were described in Appendix 
A (see the equations (A2)-(A4)). Remember that e is a small parameter of the 
perturbation theory. The scalar field operator \ which describes particles with 
mass m x can be decomposed into spatial Fourier modes 

so that every mode corresponds to a quantum oscillator with Hamiltonian (Bl), 
complex coordinate Q = Qk and frequency Q = Qk given by 

n 2 k = ul + eg(urt), (B17) 

where u 2 . = rnt + k 2 . Because the frequency Qk depends on time, the oscillator 
with label k will be excited and this means particle production. We can use all the 
above expressions obtained for the generic oscillator by simply adding a subscript 
k. Thus, the theory of scalar particle production in the case under consideration is 
only a slight modification of the oscillator formalism described above. One only has 
to take into account correctly the presence of infinite number of oscillators labelled 
by wavenumbers k and interrelations between them. The expression analogous to 
(Bll) relating the initial scalar field vacuum to the vacuum at the moment of time 
t will look like follows 

|0o>^n 7 =- P (f|4W^w)|0«>, (B18) 

where ol(i) and %(£) are the time-dependent creation and annihilation operators 
which diagonalize the scalar field Hamiltonian at the moment t. From the expres- 
sion (B18) it can be seen that scalar particles are created in pairs with opposite 
wavenumbers k. The average number of particles produced at the moment t in the 
fc-th mode will be given by the expression similar to (B12) with the only modi- 
fication being the addition of subscript k. To determine the average number Nj- 
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of particles produced we can use the formulas (B12) and (B15). We only need to 
know the solution for Q\ {t). This function obeys the equation 

Q { - ] + (ul + eg(vt))QP = 0, (B19) 



which is just the equation (Al) which we considered in the Appendix A. Initial 



conditions for Q, stem from the initial vacuum conditions (BIO) and are 



1^(0)1 = ^=^y {^QV-^QV)\. „ = 0. (B20) 



t=o 



The solutions to the equation (B19) in the resonance band were studied in the 
Appendix A. According to (A5) we introduce the value A by 

2 



•jJJ 



(§") + ^, (B21) 



where we replaced p/q in advance by the half-integer s/2 according to the results 
of the Appendix A. Now the solution to Q k will be given by a complex linear 
combination of the solutions of type (A7), namely 

Q { ~ } = -j== ( ai cos ii 1 + ia 2 cos ii 2 ) + 0(e), (B22) 

y/2uJk 

with 



^1,2 = 1^ + 6*1,2- (B23) 



The initial conditions for 01,2 and Q\p, stem from (B20) 



ai = 02 = 1, 02-01 = |- (B24) 

.1) 



If in analogy to (A15) we introduce the variables x\^ and y\^ as 

X\ = O1COS01, X2 = O2COS02, 

y\ = a\ sin 9\ , Vi = «2 sin 2 , 

then to first order in the small parameter e the equations for x\ 2 and y\ 2 will be 



(B25) 



1) The constant value of the phase a n defined in (A4) can be eliminated without loss of 
generality by a shift of the values of 61,2- 
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identical to (A16). The initial conditions for these functions follow from (B24) 

zi(0) = j/ 2 (0), x 2 (0) = -|/i(0), 

(B26) 

i 1 2, i I ' i iZii 2 1 

Fil + |yi| = F2I + |j/2i = i- 

To the lowest order in e, the expression for the value /?& is given by 

Pk = ^(xi + ix 2 + i(yi + iy2)) ■ (B27) 



It obeys the equation 



h = V+Pk, (B28) 



which follows from the system (A16), with /i + given by (A17). The initial condi- 
tions for the value /3& stem from (B26) and are 

/5*(0) = |/3 fc (0)| = ^+|a=o- (B29) 

Solving the linear system (B28) for the function (3^ with the initial conditions (B29) 
we obtain the expression for the mean number of the particles produced 

Nk = \Pk\ 2 = l _tf/\ g \2 smh 2 (/i+t) , (B30) 

where as we recall, the value //+ is given by (A17). 

The analysis developed in this Appendix can be also used in the case when the 
value of /i + is not constant in time but changes adiabatically. Namely, if 

l/i+l « /4, (B31) 

then the mean number of the particles created can be approximated by 

N k ~ sinh 2 ( / /2 + dt) , (B32) 



where the integral in the argument of the hyperbolic sine is taken over the time 
during which resonance condition (A18) is satisfied for the mode with wavenum- 
ber k, and the prefactor similar to that of (B30) has been set to unity which is 
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reasonable because at the edge of the resonance waveband (when A m \g s \) the 
adiabaticity condition (B31) ceases to be valid. 
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